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We introduce a simple model of heterogeneous catalysis on a disordered surface which consists of two 
types of randomly distributed sites with different adsorption rates. Disorder can create a reactive 
steady state in situations where the same model on a homogeneous surface exhibits trivial kinetics 
with no steady state. A rich variety of kinetic behaviors occur for the adsorbate concentrations and 
catalytic reaction rate as a function of model parameters. 
PACS numbers: 05.40.+j, 68.10.Jy 82.20.Mj 



In modeling heterogeneous catalysis, it is often as- 
sumed that the catalytic surface is spatially homogeneous 
and that all surface sites have the same adsorption en- 
ergy [1]. However, real surfaces exhibit some degree of 
spatial heterogeneity [2] which may be attributed to sur- 
face imperfections (dislocations, defects, . . .), subsurface 
effects which modify interactions at the surface, or the 
arrangement of adsorbates on the surface and induced 
adsorbate-surface interactions. 

In this Letter, we consider the effect of quenched sur- 
face imperfections on the reaction kinetics of the idealized 
monomer- monomer model [3]. Here A and B particles 
impinge upon a surface, with respective rates k^ and ks, 
and adsorb onto vacant sites. Nearest-neighbor adsorbed 
AB pairs react and desorb with rate k r , leaving behind 
two vacancies. When all surface sites have the same ad- 
sorption characteristics, this model has simple kinetics 
which can be obtained exactly by a mapping onto a ki- 
netic Ising model with mixed zero-temperature Glauber 
and infinite-temperature Kawasaki dynamics [4]. For 
\tA ^ ks, the adsorption imbalance leads to the surface 
quickly becoming saturated by the majority species and 
the reaction stops. For k^ = ks, there is a fluctuation- 
induced coarsening of the surface into growing A and B 
adsorbed islands. The concentration of reactive interface, 
or nearest-neighbor adsorbed AB pairs, Cab^), vanishes 
with time as t -1 / 2 for a one-dimensional "surface" [4] 
and as (lnt) -1 in two dimensions [5,6]. When islands 
become of the order of the size of the surface, there is 
no longer available AB reactive interface and the reac- 
tion again stops. This saturation time is proportional to 
L 2 for an _L-site chain and proportional to L 2 In L for an 
L x L surface. Thus there is no reactive steady state on 
a homogeneous surface. 

Our primary result is that surface disorder can create 
a reactive steady state. This state is different in char- 
acter from that which arises in the monomer-dimer or 
ZGB model [7] for CO oxidation on metal surfaces. Fur- 
ther, a rich variety of kinetic behavior occurs which stems 
from the interplay between the surface disorder and the 
relative adsorption rates of the two species. Our model 
can be schematically represented by the following kinetic 



steps: 



A + S^A. 

B + s"-^4b. 



As+Bg^ABt+ZS. 

The new feature is quenched surface disorder so that the 
adsorption rate of each species depends on the surface 
site, as well as on the incident particle fluxes (Fig. 1). 
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FIG. 1. Schematic illustration of the elemental processes 
on the disordered surface. Here + sites are represented by 
round indentations which preferentially adsorb (round) A par- 
ticles, while — sites (square indentations) preferentially ad- 
sorb (square) B particles. 

We choose dichotomous disorder for the surface sites: 
one type, +, with concentration c+, prefers to adsorb 
A, while the complementary — sites, with concentration 
c_ = 1 — c_|_, prefer to adsorb B. If the relative fluxes 
of A and B are p and q = 1 — p, then the disorder is 
defined by the rule that a vacant + site becomes occu- 
pied by an A with probability p + e and by a B with 
probability q — e, whereas a vacant — site becomes oc- 
cupied by an A with probability p — e and by a B with 
probability q + e. The strength of the adsorptive dis- 
order is quantified by the parameter e which obeys the 
constraint < e < min(p, q). For simplicity, we con- 
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sider the reaction-limited process with adsorption rates 
much larger than the reaction rate (kA t s , ks,s ^> k r ). 
Thus whenever a nearest-neighbor AB adsorbate pair 
reacts and desorbs, the two vacancies are immediately 
refilled. We have checked that similar behavior arises in 
the adsorption-limited process (ltA t s, &b,s <C k r ) and an- 
ticipate that our qualitative conclusions should still be 
valid for arbitrary reaction and adsorption rates. 

To develop an understanding of the kinetics, consider 
first the mean-field limit. Define x = n\/N and y = 
n^/N as the respective concentrations of adsorbed A on 
the + and — sites, where N is the total number of surface 
sites. The respective concentrations of B on the + and 
— sites are then c + — x and c_ — y. The corresponding 
rate equations are, 

x = (c+ - x)(x + y)(p + e) - x(l - x - y)(q - e), 

y= (c_ - y)(x + y){p - e) - y(l - x - y)(q + e). (2) 

Here the reaction rate k r is absorbed into the time vari- 
able and overdot denotes the time derivative. In the first 
equation, for example, the first term accounts for the gain 
in x due to a B on a + site (which is part of a nearest- 
neighbor AB pair with the A on either a + or — site) 
being replaced by an A. 

For the homogeneous surface (e = 0), the surface 
quickly becomes filled by the majority species when 
p ^ 1/2, while for p = 1/2, the A and B concentrations 
remain static. For the disordered surface there are three 
steady state solutions: two trivial, x = y = and x = c+, 
y = c_ , corresponding to B- and A-saturation, and a 
nontrivial solution, x = x co {p,e,c + ), y = y co {p,e,c + ), 
corresponding to a reactive steady state [8]. Only one 
of these solutions is stable in the long-time limit as a 
function of the parameters (p, e, c+). Their relative sta- 
bility is conveniently visualized within two-dimensional 
subspaces of the full (p, e, c+) phase space (Fig. 2(a)). In 
Fig. 2, the regions of A-saturated, 5-saturated, and reac- 
tive states are denoted by A, B, and R, respectively. The 
phase boundaries are generally described by cumbersome 
expressions [8], but in the symmetric case c + = 1/2, the 
relatively simple formula holds 



q(q ~ P \ l -<p< 1 -, 
2 ' 3~2' 



(3) 



for the B-R boundary (with a symmetric expression for 
the A-R boundary). For the temporal behavior, it is eas- 
ily found that the system approaches either saturation 
or a reactive steady state exponentially in time, except 
along the phase boundaries. However, along the A-R 
and B-R boundaries, the appropriate density asymptot- 
ically decays to zero as t _1 [8]. 



0.5 




0.5 

(a) 



0.5 



A b 


a \ 


/ B \ 


/ A \ 



0.5 

(b) 



0.5 








0.5 



1 



(c) 

FIG. 2. Phase diagram of the disordered 
monomer-monomer model: (a) the mean-field approximation 
when c+ = i (solid lines) and c+ = j (dashed lines), (b) 
one dimension, and (c) two dimensions. The phases are indi- 
cated by bold letters. Exact phase boundaries are indicated 
by solid lines, while numerically determined phase boundaries 
are indicated by dashed lines. 



For the disordered one-dimensional surface, several 
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new regimes of behavior arise whose qualitative char- 
acters are independent of the values of c + and c_ 
(Fig. 2(b)). For p + e < 1/2, B is preferentially ad- 
sorbed on both the + and — sites and the concentration 
of adsorbed A, CU(t), decays exponentially in time. Sim- 
ple bounds can be given for this decay. A lower bound 
arises by replacing the disordered system by a homoge- 
neous one with A adsorption probability equal to p — e 
at all sites. Similarly, an upper bound arises by consider- 
ing a homogeneous system with A adsorption probability 
equal to p + e < 1/2. Since C^(t) relaxes exponentially 
to zero with time in these two models, we conclude that 
p + e < 1/2 indeed defines the B-phase of the disor- 
dered system with C^(t) ~ Cab^) ~ exp(— at) with 
a = a(p,e,c+). 

For p = e = 1/4, — sites are permanently occupied 
by B, while the A and B adsorption probabilities both 
equal 1/2 on the + sites. The system therefore decou- 
ples into a succession of independent chains, each with 
fixed B boundaries. The probability distribution for the 
length L of these chains is P(L) ~ exp(— L/Lo) with 
Lo = — l/ln(c + ). From the diffusive nature of the cat- 
alytic reaction [4] , each of these chains will saturate with 
B exponentially in time with a characteristic time of the 
order of L 2 . Upon averaging over the distribution of 
chain lengths, the average concentration of A is there- 
fore given by 

C A(t) ~ Jo°° eM-L/Lo) exp(-Dt/L 2 )dL, 

~ exp (-const, x (i/ig) 1 / 3 ) . (4) 

This decoupled chain problem is also equivalent to a one- 
dimensional random field Ising model with external field 
hi = oo on a fraction c_ of the sites and hi = on the 
remaining c + = 1 — c_ sites. For Glauber dynamics, 
this problem has been solved by Forgacs et al. [9]. In 
principle, their approach could be extended to the mixed 
Glauber and Kawasaki dynamics of our catalysis model 
and would permit exact evaluation of the coefficient of 
the above stretched exponential. 

On the B-R boundary e = 1/2 — p (with 1/4 < p < 
1/2), the + sites still have a probability of 1/2 to ad- 
sorb both species, while the — sites have a probability 
< 2p — 1/2 < 1/2 to adsorb A. Because of the global 
adsorption bias, the system still saturates to B. The 
lower bound C^(t) oc exp(— at 1 ! 3 ) can b e immediately 
constructed by considering a system where A does not 
adsorb on — sites. While the construction of an upper 
bound is non-trivial, numerical simulations on a chain of 
10 6 sites show that the decay of C^(t) on this line still 
has the same form (4). Thus the B-R boundary belongs 
to the B-phase, although the asymptotic decay exhibits 
a stretched exponential behavior. 

On the line l/4<e = p<l/2, the — sites are still 
frozen with B, but there is a bias toward A on the + sites. 
This system is again equivalent to an ensemble of de- 



coupled Ising chains (with mixed Glauber and Kawasaki 
dynamics), each with — spins on the boundaries in a 
positive magnetic field. It is clear physically that the 
characteristic saturation time tl for a chain of length L 
grows as exp(_L). A Lifshitz argument of the type given 
in Eq. (4) then suggests a non-universal power law decay 
for C'A(t). A rough estimate for the behavior of tl can 
be obtained by approximating the dynamics of a single 
chain with that of a biased random walk in the interval 
< x < L with reflection at x = and absorption at 
x = L, and with bias toward the reflecting barrier. The 
survival probability of this random walker is asymptoti- 
cally given by exp(— t/ri,), with tl = (D/v 2 ) exp(vL/D). 
To make contact with the disordered catalysis model, the 
bias should be given by v = Ap — 1 and the diffusion co- 
efficient should be D = Ap(l — 2p). With this kinetics 
for a single chain, the Lifshitz argument for the aver- 
age over all chains gives C^(t) ~ t -a (P> c +) x \nt, with 
a(p,c + ) = ln(c+)2p(l — 2p) / (I — Ap). This naive ap- 
proximation agrees well with our numerical simulations 
[8]. 

For 1/4 < p < 1/2, 1/2— p < e < p (region b in the 
phase diagram) , the average adsorption rate globally fa- 
vors B but + sites adsorb A with a probability greater 
than 1/2. This dichotomy again leads to a power law de- 
cay of CA(t) ~ t~P , but with a non-universal exponent 
/3(p, e, c_|_). By comparison to the results for the case 
p = e, where the — sites are permanently occupied by 
B, the bound /3(p, e, c+) < a((p+ e)/2, c+) can be estab- 
lished. We have numerically determined the exponent [3 
for various parameter values, and a qualitatively similar 
parameter dependence as the exponent a is found. We 
define this region of the phase diagram as the b-phase to 
indicate that it still corresponds to B saturation but with 
power-law rather than exponential kinetics (as in the B 
phase). 

Finally, on the line p = 1/2, the A and B concen- 
trations are marginally stable and a saturated state is 
reached by the formation of a coarsening mosaic of ap- 
proximately equal size A and B domains. (Here we as- 
sume c_|_ = c_ = 1/2, otherwise the a-b boundary will 
be curved.) However, the disorder drastically slows 
down the interface evolution compared to the homoge- 
neous system. To appreciate the effect, consider the mo- 
tion of a single AB interface. In a region with contigu- 
ous + sites, the interface experiences a bias to the right, 
while in a contiguous region of — sites, the interface is 
biased to the left. Thus there are effective potential wells 
for AB interface motion with minima at points where a 
H — nearest-neighbor pair of sites occurs. Thus the in- 
terface motion is analogous to Sinai diffusion [10]. and 
the density of interfaces should decay as However, there 
are microscopic differences in details between Sinai diffu- 
sion and our disordered catalysis model, and the interface 
density is numerically found to decay as 
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C A B(t)~(lnt)-i^, (5) 

with a non-universal exponent 7(e), in contrast to 7 = 2, 
if the analogy with Sinai diffusion was exact. 

For a two-dimensional disordered surface, five phases 
occur, including a reactive steady state R. As in one 
dimension, it may established rigorously that in the A 
phase, where A is preferentially adsorbed on both the 
+ and — sites, A-saturation is approached exponentially 
in time. This conclusion is independent of the values of 
c_|_ and c_ . In numerical simulations on a square lattice 
of 10 6 sites, we have also found a and b phases for the 
symmetric case c + = c_ = 1/2 (Fig. 2(c)). In these 
regions the minority species decays as a power law in 
time, C m i n ority(^) ~ t -13 " 1 , with a nonuniversal exponent 
f3 m = /3 m (p, e, c_|_). Finally, in the reactive steady state 
(R) the densities approach nonzero values exponentially 
in time. 

More subtle kinetics arises on the phase bound- 
aries. On the a-A and b-B boundaries the minority 
species appears to exhibit a stretched exponential decay, 
ln(C m i n ority(^)) ~ —t Sm - To obtain a theoretically plausi- 
ble estimate for this exponent, it is convenient to consider 
the extremal point p = e = 1/4 (or q = e = 1 /4) and fol- 
low the same line of reasoning as in one dimension. For 
p = e = 1/4, for example, the — sites are permanently oc- 
cupied by B, while on the + sites, the A and B adsorption 
probabilities both equal 1/2. The — sites thus play the 
role of "sources" of the 5-phase and the long-time behav- 
ior is determined by the filling of large voids which do not 
contain — sites. Within the voids the kinetics is just that 
of the homogeneous and unbiased monomer-monomer 
model (i.e., p eff = p + e = 1/2, q eff = q - e = 1/2). 
By again appealing to the mapping onto onto a kinetic 
spin model [5] , it may be argued that the average concen- 
tration of A in a void decreases as exp(— t/R 2 ). Since the 
probability distribution for disk-shaped voids of radii R 
is Poissonian, P(R) ~ exp(— R 2 /R 2 ), the configuration 
average A concentration is 

C A (t) ~ exp - ±^dR ~ exp (-W 2 ) . (6) 

Therefore S m = 1/2 at the extremal points. Based on the 
equivalence between catalysis kinetics and the trapping 
problem in two dimensions, we expect the universal value 
S m = 1/2 along the a-A and b-B boundaries. 

Finally, on the a-R and b-R boundaries we observe 
very slow decay of the minority species. It is difficult to 
establish whether this decay is a power-law with tiny ex- 
ponents or logarithmic. The logarithmic behavior seems 
more natural and qualitatively agrees with the (lnt) -1 
decay for the number of interfaces which is expected to 
occur at the extremal point p=l/2,e = (see [5]). 

In summary, we have found a reactive steady state in 
the monomer-monomer model of heterogeneous catalysis 



on surfaces with adsorptive disorder. This result sup- 
ports intuitively appealing hypothesis that disorder gen- 
erally stabilizes the reactive steady state. We have also 
observed rich kinetics for the reactant concentrations and 
reactivity in our catalysis model, including exponential, 
stretched exponential, power law, and logarithmic decays 
with time for different system parameters. 
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